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ABSTRACT 


Model tests at M.I.T. have indicated that the cam- 
type wave power generator is capable of high efficiencies 
of energy extraction from ocean waves. This report presents 
a design methodology for determining costs of energy produced 
from wave power generators. Feasible deSigns were developed 
for ocean locations around the coast of the United States 
with the resulting breakeven costs with oil ranging from 
14 to 30 S/BBL. 


Thesis Supervisor: A. Douglas Carmichael 
Title: Professor of Power Engineering 


Thesis Advisor: Stephen C. Graves 
Title: Assistant Professor of Management Science 








Siete ee 


a 


ACKNOWLEDGEMENTS 


I would like to express my gratitude to Professor A. 
Douglas Carmichael who first kindled my interest in the 
subject. His thoughtful advice was most appreciated. 

Thanks also are due to Professor Stephen C. Graves, 
Mr. Jim Mays and Mr. Noel "Brye" Davis for their helpful 
suggestions, and to Mrs. Sandy Margeson who typed the 


Manuscript. 





ABSTRACT. 


TABLE OF CONTENTS 


PeRNOWLEDGEMENTS. . « © « © «© © © © © © © © © © © « 


Meee OF CONTENTS . . -. « «© © © © © © © © © © © «© «© 


foo T OF SYMBOLS e e e e e e @ e e e e e e e ® e e e 


Chapter l1: 
Chapter 2: 
Chapter 3: 
Bie : 
BZ: 
a : 
Chapter 4: 
Gaapcer 5: 
Chapter 6: 


Gel 


OCMC ETOM. Ueueliss 6 6 « «© «© 0 « « « 
Theoretical and Experimental Background 
PME WAVES 5 4 6 6 + «© +6 « « 6 @ « 
MNyeeorsecrical DeErilVation. s+ . .« « « » « -« 
Limitations on Wave Power ....... 
Prepabilidity wPaSeFIDUCLONS . 1... « « « + 
WaviemelmIMatOlOgy. .< =. « « <« « « © w« « « 
Formulation of Model for WPG Evaluation 
Simulation Results and Analysis .... 
Cam Diameter - No Separation. .... . 


Minimum Cost/KW Diameter Selection - No 
Separation. e® e e e e e e e e e e@ e@ ® ® 


Variability of Wave Period. .... . 


CanlseparatVOn. =o. . see ew lt lt ls 


Minimum Cost/KW Diameter Selection - With 


Bielelinchel@ Qi ute. « + » « 6 © «© 6 5 6 « 


Energy Storage, Power Rating and Availability 


Optimal Level of Energy Storage .... 


Variability of Wave Power ...... . 


Bu 
52 
oe 


62 
68 
63 
qe 





@hapter 7: Feasibility Tradeoff Studies. 


7.1: Seasonality Effect on 
Power Level ee at xfs 


Pepe Structure ...... 
Pees Conversion. .... - 
7.4: Electrical Generation 
7.5: Optimal Power Rating. 
feo: Transmission. .... 
eee Energy Storage. . . 
7.8: Other Subsystems... 
Chapter 8: Cost Analysis .... 
memes, §6COStTS «2 Ow te tl ltl lt 
See: Power Analysis. .. . 
Meee Availability. ... . 


8.4: System Design... . 


Selection 


8.5: Design Method for Fuel Saver Analysis . 


8.6: Design Method for Baseload Analysis. . 


8.7: Hypothetical Ocean Energy Farm Cost 


Mameulation ... .« .« 


8.8: Value of Energy Storage... 


See Structural Cost Reductions. .. 


SeeeuseAmalvysis of Costs . . 


Chapter 9: Conclusions and Recommendations 


MP RRENCES. . . . « «© « « 
Appendix A: Model Description 


Appendix B: Program Listing 


14 


74 
73 
oil 
34 
86 
vie, 
100 
102 
106 
106 
OT 
110 
ia 
teas: 
bilge 


oe 
2s 
La3 
133 
32 
139 





a Q Wt VU KR TF 


+3| 


PoreOr wor MbOLS 


wave amplitude 

wave phase velocity 

wave group velocity 
deepwater group velocity 


2A - wave height 


= 72H ~ significant wave height 
S rms 


cam diameter 

cam diameter at maximum power out 
wave energy 

gravitational constant 

water depth 

wave number 

wave power 

mean wave power 

effective length multiplier with separation 
wave period 

mean wave period 

Separation distance 

water density 

wavelength 


frequency 





observed wave height/mean wave height 
standard deviation of wave period 
standard deviation of wave power 


efficiency 





CHAPTER 1 


INTRODUCTION 


Energy from renewable resources offers an alternative to 
conventional fuel-fired energy production systems. Escalating 
fuel prices, political/economic desires for energy indepen- 
dence and reduction in environmental pollution are all 
incentives for developing renewable resources. The cam-type 
wave power generator is one such resource, utilizing the 
energy in ocean waves. 

The idea of extracting the energy in ocean waves by a 
cam-Shaped nodding duck was first presented by Stephen Salter 
of the University of Rainburgh iel974 (ret. 26). Since 
waves drive each other with high efficiency, Salter believed 
that an optimally shaped device would cause an incoming wave 
to react as though the device was just another wave. 

In deepwater travelling waves each particle moves in a 
Meably circular orbit. At the surface the diameter of the 
orbiting particles equals the wave height. Below the surface 
the orbital diameters fall off exponentially with depth. 

Salter's duck shape is displayed in figure 1-l. The 
wave side of the cam was designed to conform as nearly as 
possible to the diameters of the orbiting particles. The lee 
Side of the cam was designed as a cylinder because this shape 
has constant displacement during the rocking of the cam and 


therefore transmits no energy. 








meecune tls Profile of Salter Cam. 
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Salter tested a model of his device and found that it 
was capable of extracting nearly 90% of the wave energy inci- 
dent on its length. Experiments conducted at M.I.T. with a 
similarly shaped cam also yielded high efficiencies. 

Wave energy, like all other renewable resources, is 
characterized by low intensity and random flux. In order to 
predict the useful absorption of a wave energy device, data 
must be collected and analyzed for the desired ocean location. 
The following parameters adequately describe ocean waves for 
purposes of energy collection: distributions and mean values. 
of wave height, wave period, variability of wave period, power 
density and variability of power density. 

Before its economic feasibility may be determined, the 
Wave Power Generator (WPG) must be measured for performance 
in random seas. To accomplish this a model must be developed 
which is capable of dealing with environmental variables like 
wave height, wave period, period variability, power density 
and power density variability, and with design variables like 
cam diameter, energy storage, power rating and separation 
between WPG units. 

Once performance is determined, subsystem components 
for structure, conversion, electrical generation, transmission 
and energy storage may be selected on the basis of feasibility, 


cost and efficiency. 


iil 


The estimates of deliverable power and subsystem costs 
enable determination of the economic viability of the WPG 
with respect to conventional power sources. 

This report deals with the four basic steps outlined 
above: 

(1) Collection and analysis of ocean wave statistics; 

(2) System performance modelling; 

(3) Subsystem component selection; 


(4) Analysis of investment opportunities. 


=— 
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CHAPTER 2 


THEORETICAL AND EXPERIMENTAL BACKGROUND 


Experiments on scale models of the cam have been con- 
ducted at M.I.T. and by Stephen Salter of Scotland. The basic 
characteristics of the M.I.T. cam are depicted in figure 2-l. 
This cam was tested at M.I.T. in both a 2-D mode and a 3-D 
mode: in the 2-D mode the cam's length filled the width of a 
narrow tow tank so that waves generated in the tank had a 
crest length equal to the cam's length; in the 3-D mode the 
cam was tested in a larger tow tank allowing separation 
between the ends of the cam and the tank's walls. In both 
cases the cam was rigidly fixed in space, being allowed only 
to roll about the axis of its length. This axis is parallel 
meewave crest. 

The experiments at M.I.T. were conducted by N.B. Davis 
under the direction of Professor A.D. Carmichael and the 
results are documented in a thesis by Davis (ref. 7). 

“The efficiency of the cam as a wave power absorber is 
defined to be the percent of power in a wave of crest length 
equal to cam length which the cam absorbs. Davis tested the 
cam's efficiency as a function of non-dimensional frequency, 
wia/g, for three levels of inertia: low = .060 Kam“, middle = 
2063 Rea and high = .081 Lone His resulting efficiency 


curves are graphed in figure 2-2. Also included in this 
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Figure 2-1: M.I.T. Cam Geometry 
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figure is the efficiency curve published by Stephen Salter. 
In each of these curves there is a fairly wide band of high 
efficiency levels. IncreaSing inertia appears to lower the 
frequency of optimal performance. 

When the cam was tested in the 3-D mode, certain fre- 
quency waves produced very large increases in efficiency. 
This efficiency curve along with the associated 2-D curve for 
the cam with middle inertia is plotted in figure 2-3. Besides 
increasing the cam'’s efficiency, the 3-D mode displayed 
optimal efficiencies at lower frequencies than the 2-D mode. 

Davis defined a non-dimensional cam response parameter, 
S9a/H, where 8 = peak-to-peak radian angular displacement, 

a = cam radius and H = wave height. He plotted response as a 
function of frequency for various levels of damping and this 
plot is included in figure 2-4. Superimposed on this plot 
are dashed lines of efficiency which show that the maximum 
efficiency occurs at the damping level which causes the 
response (§a/H) to have a value ~ 1.0. The level of damping 
was found to have no effect on the frequency of optimal 
performance. 

In a paper by K Budal (ref. 2) it is shown that the 
idealized Salter cam construction comprised of a linear row 
of closely spaced cams is capable of absorbing 100% of the 
power incident on its projected length. This helps to 


explain why the 2-D efficiencies in figure 2-2 are so high. 
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Budal also describes the possibility that a wave interaction 
effect on wave-absorbers separated by a distance of approx- 
imately one wavelength (which is large compared to the 
absorber's length) will cause dramatic increases in the 
efficiency of each absorber. This may explain why the M.I.T. 
experiments showed an increased efficiency when the cam was 


operating in the 3-D mode. 
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CHAPTER 3 


ENERGY IN WAVES 


B.1 Theoretical Derivation 

The basic requirement of the wave power generator is to 
meeeure the energy that is transported across the surface of 
the water in the form of waves. Given a particular set of 
wave characteristics the level of power per unit length of 
free surface (kilowatts/meter) that exists in a train of waves 
may be determined. 

Wave power is the product of the energy in a wave times 
the rate at which the energy is being transported. For plane 
progresSive waves the average energy per unit length of free 
surface is 

E = 50 ga‘ (1) 

The average rate of energy flux of a plane progressive 


wave across a vertical control surface is 

2 
gwA” /K (2) 
The rate at which the wave energy propagates is the 


mean energy flux rate (dE/dt) divided by the mean energy per 


unit length (FE): 








20 
Equation 3 shows that the energy in a wave propagates 


at the group velocity WE a < Since the power in a wave is the 


peocuct of energy times rate of energy propagation, 


af 
Te © gH, C os 


For all depths of water, 
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The value of [= + | may be graphed (figure 3-1) 
2 SinH (2Kh) 


or tabulated (table 3-1). 
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TABLE 3-1 


VALUES OF an FOR VARIOUS WATER DEPTHS AND WAVE PERIODS 


Water Depth (M) 


S 6 S) 1L5 20 30 
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B. 2 Limitations on Wave Power 

Before examining actual wave power Statistics from 
various sites it will be worthwhile to survey some of the 
environmental phenomena which limit the development of power. 

Wind speed is the most important factor in the develop- 
ment of waves. Darbyshire, Sverdrup-Munk and Neumann have all 
proposed spectra relating significant wave height to wind 
speed (U). These spectra show that He Lome renorilonal to 
ue or aoe (ref. 15). This implies that the highest values 
of wave power will occur in regions or seasons with very high 
continuous wind speeds. 

Fetch, the distance over water that the wind blows in 
an unobstructed manner, is another important element in wave 
power limitation. For a given wind speed significant wave 
height increases as fetch increases. Figure 3-2 displays 
the effect of fetch on Ho: 

Related to fetch is sheltering. Harbors and other 
sheltered waterways cannot possibly experience the continuously 
high levels of power available in the open ocean. 

In very shallow water locations wave energy may be 
dissipated in wave-breaking and bottom losses. 

This very general survey should aid in predicting 
locations of the highest power densities. Since winds pre- 


vail from the west for the U.S. it is reasonable to expect 
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coastal power densities to be higher on the West Coast than 

on the East Coast, because of the effects of fetch and 
sheltering. Higher power densities are anticipated at deep 
ocean locations rather than coastal locations because of the 
effects of fetch, sheltering and dissipation. Since winds 

are more severe in the winter months higher power densities 
are expected during winter. The actual values of wave power 
at different locations and during different seasons are 
detailed in Chapter 4 and they substantiate these very general 


observations. 


SAS) Probability Distributions 
3.3.1 Wave height. In spectral analysis of wave 


records the Rayleigh distribution is generally assumed for 
wave heights (refs. 4, 15). The simplified form of this 
distribution is 


2 
ope Se pera (11) 


where € = observed wave height/mean wave height. 

A statistical analysis of this distribution was per- 
formed on 1963 samples from statistics in ref. 12. The 
Rayleigh expression was inverted using the computer library 


program MIT-SNAP and a regression line was fitted to the 





ay 


inverted distribution. The regression yielded + = oo 7e wile 


small and evenly distributed residuals. Although the ot is 
not particularly high, the regression does indicate that the 
Rayleigh assumption is a reasonable approximation to the 
observed data. 

3.3.2 Wave period. If little is known about the wave 
period at a particular location then a good approximation 


(ref. 10) to the distribution can be achieved by using 


T, 4 
a7 675(=) 


P(=) =1- = 


ly 


(12) 


This equation yields a distribution which is symmetric about 
its mean. 

The statistics used in estimating power in this report 
all assume that periods are normally distributed with a 
standard deviation On A regression line was fitted to the 
inverted normal distribution for the same 1963 samples used 
in the wave height analysis. The value of r? was .96 and 


the residuals were small and uniformly distributed. 


Seoes Independence of height and period. Wave heights 


and periods are generally assumed to be independent (refs. 20, 
32). A statistical test for correlation between height and 


period was performed using 1000 observations from joint 
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height-period tables scattered over all seasons at Nags Head, 
North Carolina (data from ref. 32). Mean wave periods for 
each of ten wave height groups were tabulated an a hypothesis 
test for correlation was conducted using a 5% risk of type I 
error. Since the data was reduced to only ten wave height 


groups the test was performed using manual methods and the 


Serstribution. 


H, =: 0 x 0 





Test for t = ya 


1-r* 
Bere r* = explained variance _ i - s* (T/H) 
total variance — 52 
a 


n = number of observations 
The assumed form of the regression equation is 
par/H) = & + BH 


s°(t/H) = 2(rr? - @ rT - BITE) 
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s*(T/H) = 77(822.58 - 8.81 (90.4) - .041 (500.6)) 
= 868 
eee 555 9404) 00.4)) = .536 
tT = I0 
r* = E aS .05 
r= .224 
t= .224 /-or, .65 


For a 5% risk of Type I error, t = 2.3. Therefore, accept 


hypothesis H 0 = 0, that there is no correlation between 


0 
mand T. 


3.3.4 Wave power. Nath (ref. 20) derived the pro- 
bability distribution for wave power under the assumptions 
that Power (P) ~ Hq“, which has already been shown to be 
valid, that height can be represented by the Rayleigh distri- 
bution and that i can be represented by the Rayleigh 
distribution. This last assumption for p? is the same 


distribution presented in section 3.3.2 for use when wave 


@eriOd Standard deviation is unknown, that is, 
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T, 4 
pi) 2 - 9 ~ 67512) 
al 
2 
Let y = H T = power 
y = mean power 


Using the above assumptions, Nath derived the 
following distribution: 
-y/y 


P,\Y) = e 


L 
y 
y = standard deviation 


A statistical test of this distribution was performed 
using the same 1963 samples that were used in the height and 
period tests. The distribution was inverted and a regression 
line was fitted to the inversion with 1? = .99. The regres- 
sion strongly suggests validity of this exponential distribution. 

A further test was conducted on standard deviations of 
wave power arrivals. The exponential distribution yields an 
Me pected value of 0,/P = 1.0. Data from Pierson (ref. 22) 
was utilized with 12 monthly summaries at 20 different loca- 
tions. A regression line was fitted on the observed standard 
deviation of power (01) Vsmeene cieoretical standard 


deviation (P). The regression line was 





SZ 


d, = 2.83 + 1.03 P 
with x? as O06 

With values of Op ranging up to 350 KW/M in the data, 
the regression constant (2.83) is negligible. The slope of 
1.03 compares favorably with the expected slope of 1.0. 
These results further validate the use of the exponential 
power arrival distribution. 

It is important to note that the time between observa- 
tions of wave power may play a significant role in the shape 
of the distribution. Most of the records which were used in 
this report are based on monthly summaries of observations 
taken every 12 hours. The Ace nueion of power over any 
given 12 hour interval is clearly not as severe as the 
exponential distribution, largely because of the persistent 
nature of weather patterns. 

It is beyond the scope of this report to attempt to 
determine the shape of the power distribution over a 12-hour 
interval, but a confidence interval about the expected value 
of 0,/P can easily be established. Using 150 sets of data 
from the Marine Environmental Service, Canada (ref. 38), with 
each data set containing 5 observations of wave energy in a 
given 12-hour interval, a mean and standard deviation of 


energy were computed. 
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n = 150 

x = sample mean = O/E =a h56 

S = sample standard deviation = .132 
eet = Ce =) 0108 


I+ 


x + 1.96 te .156 1.96 (.0108) 


= ees Oe. O 2h 
mamge: «135 > .17/7 


In other words over any given 12 hour interval the 
standard deviation of wave energy is expected to be 15.6% 
2.1% of the mean wave energy. This is obviously much less 


than the expected 100% using the exponential distribution. 
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CHAPTER 4 


WAVE CLIMATOLOGY 


The problem of defining ocean waves statistically is 
not new, yet, it is only within the past decade that acceptable 
techniques for recording data have been developed. The U.S. 
@emecOorps Of Engineers Civil Engineering Research Center 
(CERC) has compiled the most accurate data available by its 
system of wave gauges around the coast of the U.S. (ref. 32). 
Unfortunately, all of the wave gauges are located within one 
mile of the coast and are, therefore, not very useful in pre- 
dicting open-ocean power densities. The manual Ocean Wave 
Statistics (ref. 12) provides an acceptable compilation of 
joint wave height and period distributions for very large 
open-ocean areas around the world. These two sources and 
references 22 and 38 provided the data for analysis of ocean 
waves around the U.S. 

Since this report deals with a feasibility study, the 
statistics of ocean waves will be most useful if they present 
a very broad picture. Raw data from 17 coastal locations and 
29 open-ocean locations was analyzed and summarized to present 
as broad a picture as possible, but still to retain enough 
identity to perform calculations necessary for the study. 


Tables 4-1 to 4-5 contain the following variables: 





A, 


Region - North, Middle and South Atlantic and 
North, Middle and South Pacific off the coast 
oe ene Uso. 

Ocean location - Coastal waters (within 1 mile) 

Open-ocean (about 100 miles off shore) 

Seasons - Dec-Feb; Mar-May; Jun-Aug; Sep-Nov 

Wave Height - H in meters 

Wave Period - T in seconds 

Wave power density - P in KW/M 

Non-dimensional power variability - o,/P 


Non-dimensional period variability - O/T 


It is worthwhile to point out some of the strengths and 
limitations of these tables. Hundreds of thousands of obser- 
vations from over 300 records were summarized to create these 
tables; therefore, the tables should not be applied to any 
specific location. Data for coastal regions is more accurate 
than data for deep ocean locations. Calculations of power 
density were performed by the methods outlined in Chapter 3. 
The values of wave period (T) in the deep-ocean sections for 
Atlantic regions are probably not very accurate, as these 
Observations are based on shipborne wave-recorders. For this 
reason and since wave period remains essentially unchanged 
With depth (ref. 32) only one value of O/T was calculated 


for each region, based on coastal readings. 





The graphs in fig. 
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power level over the seasons for the East and West coasts. 
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TABLE 4-4 
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Figure 4-1: Variation over year of mean power non- 
dimensionalized to mean monthly average. 
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CHARTER GS 


FORMULATION OF MODEL FOR WPG EVALUATION 


In the previous chapters experimental evidence which 
substantiates the viability of the cam-type wave power 
absorber and the environmental phenomena which affect the 
cam's operation were presented. The next step in the feasi- 
bility study is to analyze the cam with all its associated 
parameters to determine how it might operate in a real-world 
environment. This task may be viewed in the pattern of the 


flowchart below: 








Formulation of 
the problem 

















Data & Statistics: 
Collection & 
Processing 


Techniques: 
Math Programming 
Simulation, etc. 








Model: 
~ Formulation, 
Validation and 

Experimental Design 















mormulation of the Problem 

The cam system may be characterized by length, dia- 
meter, separation between elements, natural frequency of 
oscillation, motion damping, energy storage, a scheme of con- 
version to energy in a useful form and the users' demand for 


that energy. 
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The cam's environment (the ocean) is essentially 
characterized at any location by distributions of wave 
height, period or frequency of waves, power one SASS) 10 
waves, water depth and distance from land. 

The coupling between the cam and its environment is 
characterized by sets of efficiency curves which are functions 
of cam design and the environment. 

The basic problem is the same for all energy-related 
systems: that is, to optimally design the parameters of the 
cam (at the feasibility-study level of detail) and to locate 
the system so that it delivers the least expensive energy to 
the users. The model must deal with the stochastic nature 
of the environment and be flexible enough to determine the 


effect on performance of variations in system design. 


@eeeection and processing of real world data and statistics 


Cam efficiency characteristics and environmental 


variables have been presented in previous chapters. 


Formulation of the model 

There 1S no necessity to deal with both performance and 
costs simultaneously. Performance will be analyzed first. 

The model must be capable of handling discrete-time 
random inputs, a large number of variables and parameters, 


mathematically intractable efficiency distributions, 
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comparisons of alternative systems or operating conditions, 
reproducible experiments, and output statistics in the form 
of summaries with frequency distributions. | 

The obvious model type selection for handling the 
above conditions is simulation. The simulation language GPSS 
was chosen for the following reasons: (1) the author's 
familiarity with its use; (2) discrete-time capabilities; 
(3) ability to deal with complicated distributions; (4) many 
commands for assembling statistics; (5) debugging capabilities; 
(6) transaction flow tracking capabilities; and (7) repeat- 
ability of experiments. The model description is included 


in Appendix A and the program listing in Apvendix B. 
Model Validation 


Several experiments were performed on the model to 
check its logical and statistical validity. From the table 
POAVA the distribution of power available was plotted to 
check its correlation with the input exponential distribution. 
several efficiency curves were reproduced by simulating non- 
random seas and extracting data from the table EFFIC (effici- 
ency). Statistics from the table PEROD (period) were 
correlated with several different input period distributions. 
Individual model segments were checked for proper operation 
using the block count table (a standard output with GPSS) to 
Verify that transactions (wave power arrivals) proceeded ina 
logical sequence when the model was constrained to force operation 


Of each segment. 


- _ 
Orion. 
; 
| 
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Experiments were designed and performed to test the 
effects of changing environmental conditions and system 
parameters on system output. These experiments with the 


analyses will be discussed in Chapter 6. 
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CHAPTER 6 


SIMULATION RESULTS AND ANALYSIS 


The simulation model was run with a large number of 
different configurations in order to develop design lanes or 
determine the effects on the system of variable cam diameters, 
hydrodynamic efficiency, wave period variability, wave power 
variability, energy storage, power demand, power rating, 


power availability and cam separation. 


B. 1 Cam diameter - no separation 

Iterations were performed with the model to determine 
the diameter of maximum power output for a range of wave 
period variances. The maximum output diameter is independent 
of wave period variance but is different for each of three 
levels of inertia tested at M.I.T. The results are listed in 
Table 6-1. 

TABLE 6-1 


CHARACTERISTICS OF MAXIMUM POWER OUTPUT DIAMETER 


Inertia wva/g A/D n/n 





ee 
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Note that the maximum power output occurs near the 
apex of efficiency curves from the M.I.T. experiments for 
each level of inertia. The high inertia asst delivers more 
power with a smaller cam diameter (smaller wWa/g) and there- 
fore 1s the optimal choice of the three. The efficiency 
listed in table 6-1 is non-dimenSionalized to the efficiency 
of the high inertia value. The actual number value of the 
efficiency is a function of the variability of wave period 
about its mean: this will be discussed in section 6.3. 

The question next arises as to whether or not this 
Single cam diameter indeed delivers the most power from the 
faves Since wave periods are varying about the mean and optimal 
cam diameter is a function of wave period. The model was 
tested to see if varying the cam diameter about its optimal 
mean would produce any increase in system efficiency. An 
arbitrary entity was devised representing the standard devia- 
tion of cam diameters about the mean value obtained from 


table 6-1. This entity is defined o Figure 6-1 shows the 


D* 
Simulated results of power out and standard deviation of power 
But vs. Op/ om: The results indicate that the maximum power 
Output diameter for a system of cams in a given random sea is 
a point design, in other words, all cams in the system should 
have the same diameter, since the mean power decreases and the 


Spread of power increases as diameters are increasingly varied 


about the point mean. 








Power/Max Power 


Oa. 


46 


Mean Power 
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re 6-1: Effect of Varying Cam Diameters about 
Point Maximum. 
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Since it is possible that the minimum cost/KW system 
will have a smaller diameter than the maximum output system, 
off-design performance was analyzed with the model. The 
BerLormance is depicted in figure 6-2, where power (P) is 
non-dimensionalized with the maximum power output Caer a 
and diameter (D) is non-dimensionalized with the diameter at 
which maximum power is delivered (Ds) 


The spread of power delivered changes considerably as 


diameter is reduced from Die due to the filtering character- 


x 
istics of the efficiency curve. Because of this, figure 6-2 
must be adjusted to ensure that all points have equal 
availability. The availability, or the percent of time that 
Reyer is equal to or greater than a specific power level, was 
extracted from the simulation model output tables. The 
availability of power at each point of the curve in figure 6-2 
was compared with the availability of the maximum output 

power eas = 1.0). The curve of adjustment factors is 
included in figure 6-3. 

Assuming that the marginal value of increased availa- 
bility is equal to the marginal cost of increased power, the 
adjustment factor curve may be multiplied point by point with 
the off-design power curve to yield the adjusted curve of 
figure 6-4. This adjusted curve is used when cost estimates 


are applied to find the minimum cost diameter. 
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Figure 6-2: Delivered Power vs Cam Diameter - 


No Separation 
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Figure 6-3: Availability Adjustment Factor For Off- 
Design Diameters — No Senvaration 
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Figure 6-4: Delivered Power vs Diameter Adjusted 


For Equivalent Availability. 
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or. 2 Minimum cost/KW diameter selection - no separation 
With the off-design performance of the cam known, all 
that is needed is the corresponding cost asia es find the 
minimum cost/KW design diameter. The basic trade-off involves 
only structural costs and delivered power. According to 
Kyrkos (ref. 17) structural costs are directly related to the 
weight of structure. A useful structural weight estimating 


relationship (from ref. 25) is 


W, = 1.68341 ene Gp 26701720 CN = 23.283 tons a) 
where CN = LBZ x 10> 

L = length 

B = beam 

Z = depth 


With the Salter cam, structural weight is essentially 


embodied in the circular cylindrical main body. Therefore, 


beam = depth = cam diameter and CN = LD“ x 10>. 
Bes constant length, 
fee ces4i x20 *° p- + 4671721 x 10° v2 


iL 


18, 





D2 


Therefore, off-design cam costs are proportional to p*. 
The resulting cost curve, non-dimensionalized to the cost of 
the cam at D_., is drawn in figure 6-5 along with the off- 


design power curve. Dividing the $/$ curve point-by-point 


$/$ a 
max , 
by the Be as curve yields P7P_ Or non-dimensional $/KW. 
From the figure, $/KW is minimized at D/D, wax = .60. 


In other words, the minimum cost/KW system will have a 

diameter which is 60% of the diameter at maximum power output, 
where the cam delivers 68% of the power of Das at 36% of the 
cost of Dd. 


ax 


6.3 Variability of wave period 

The non-dimensional parameter chosen to describe the 
Variability of wave period is O/T. A given cam design will 
deliver less power as the value of O/T increases, Since the 
cam wili be caused to operate outside its region of maximum 
efficiency more often with higher o,/T. The simulation model 
was run with various values of Om and various values of T 
yielding the curves in figures 6-6, 6-7, and 6-8. Figure 6-8 
represents the efficiency curve adjusted for equivalent 
ae@eaiability. 

Typical values of O/T (from the tables in Chapter 4) 


are displayed in figure 6-8 with the resulting 2-D cam 


efficiencies between 60-70%. 
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Figure 6-5: Minimum Cost Diameter Selection — No 
Separation. 
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Figure 6-6: Cam Efficiency -— No Separation. 
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Figure 6-7: Availability Adjustment Factor Curve 
For Cam Efficiency Curve. 
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Figure 6-8: Efficiency Curve Adjusted For Equivalent 


Availability - No Separation. 
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6.4 Cam separation 

Budal (ref. 2) has shown that separation may effectively 
increase the efficiency of a wave power absorber by wave 
interaction. He defines a factor, q, where 


N 


q= [1 + = (N-n)cos(nKd siny) J_(nKd)]_ 
=| ~ 


1 


where q = effective efficiency multiplier 
N = number of power absorbers 
K = wave number 
d = separation distance 


Y = angle between cam and normal to wave crest 


Davis (ref. 7) displayed the effect of separation by 
drawing efficiency curves for a 2-D mode and a 3-D HO In 
the 2-D mode maximum efficiency was about 95%, while in the 
3-D mode maximum efficiency was about 175% (figure 2-3). 
Since the 3-D mode data was taken in an open tank, the cam 
effectively "saw" images of itself caused by the two tank's 


Nee 
walls. Davis' data was reexpressed to the form g = — 5 and 


2-D 
plotted against Kd in figure 6-9. The separation distance, 





d, was taken to be twice the distance from the cam's center 


to the tank wall, or the distance from cam to image. 
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Figure 6-9: Effect of Separation Using MIT 


Experimental Data. 
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There appears to be a discrepancy between the observed 
"q" and the "q" expected from Budal's analysis. Assuming 
N = 3 and y = 0, Budal's equation yields, 


q = [1+ £(2 J.(Ka) + 3. (2Ka))]~ (3) 


This equation is plotted as the solid line in figure 
6-10. For the observed data plotted in the same figure, the 
value of Kd was divided by 3 to achieve the visible correlation 
between the observed data and Budal's equation. This is the 
discrepancy referred to previously, and it is also the reason 
that the horizontal axis in that figure is left unmarked. 
Because of the discrepancy, conclusions about the optimal 
value of Kd may not be drawn with certainty. However, it may 
be stated that there exists some value of separation distance 
which will yield appreciable increases in cam efficiency. 

With the 2-D cam, efficiency was shown to decrease 
as the value of O/T increases. This same effect occurs 
with the 3-D cam. Figure 6-11 shows the cam efficiency as 
merunction of o,/T based on a cam of middle inertia designed 


at the maximum power output, that is, at wv¥a/g = .63. 
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Figure 6-10: Separation Effect Comnarison Between 
MIT Observations And Budal Equation. 
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Figure 6-11: Cam Efficiency Curve With Sevaration 
At Optimal Sevaration Distance Adjusted For 
Equivalent Availability. 
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B.5 Minimum cost/KW diameter selection - with separation 
It was shown in section 6.2 that the optimum (minimum 
cost/KW) diameter is not the same as the maximum output 
diameter. An analysis similar to that performed on the 2-D 
cam was performed on the 3-D cam for off-design diameters. 
Davis showed that the maximum power output for the 
3-D cam with middle inertia occurs at a value of wv¥a/g which 


is lower than that for the 2-D cam: 


COMPARISON OF MAXIMUM POWER OUTPUT CHARACTERISTICS FOR 


THE 2-D AND 3-D CAMS 


wla/g A/D 
Di=)Bi 10 5.4 
IB, 563 Fad 


Table 6~2 shows that for a given frequency wave optimal 
cam diameter is smaller for the 3-D cam than the 2-D cam. The 
diameter at maximum power output is 32% smaller for the 3-D 
Sam, Or, De ee (3-D) = .68 D ax (2-D). 

Figure 6-12 shows the off-design 3-D performance; 
figure 6-13 shows the availability adjustment factor; figure 
6-14 shows the adjusted off-design performance curve and 


figure 6-15 shows the minimum $/KW diameter. The analysis is 


Similar to that performed in section 6.2. 
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Figure 6-12: Delivered Power vs Cam Diameter - 
With Separation. 
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Figure 6-13: Availability Adjustment Factor For Off- 


Design Diameters — With Sevaration. 
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Figure 6-14: Delivered Power vs Diameter Adjusted 
For Equivalent Availability - With Senaration. 
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Figure 6-15: Minimum Cost Diameter Selection - 
With Sevaration. 
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Note that the minimum cost diameter for the 3-D cam is 
mpout .52 Best while for the 2-D cam it is - 60 Denis Separa- 
tion then will play an important role in reducing costs of the 
system for these three reasons: 

1. Higher efficiency for a given O/T 

2. Smaller maximum output diameter (Dee 

3. Smaller minimum cost diameter ratio (DDE =) 


ax 

These observations may be used to provide a rough 
estimate for the value of separation. At a typical value of 
o,/T = .30, the factor reduction in diameter of maximum 
power output is .68. The factor reduction in the minimum 
cost diameter is .52/.60 = .867. Using ences three factors 


and again assuming that structural costs are proportional to 


:, the structural $/KW of the cam with separation is 


Mee) (.867)- = 
ie 


Aa) 


times the structural $/KW of the cam with separation, assuming 
no additional costs for separation. This remarkable savings 
indicates that further research should concentrate heavily in 


the area of cam separation. 
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6.6 Energy storage, power rating and availability 


A very definite trade-off exists between demand 
(assumed to equal power rating for this analysis) and 
availability (% of time that power is equal to or greater 
than demand), for any given level of energy storage. A 
family of curves of demand vs. availability may be drawn 
for any level of energy storage uSing results from simulation 
model runs. As the example in figure 6-16 displays, there is 
Mevery Sharp “knee" in each curve of demand vs. availability. 
Assuming the marginal value of increased availability equals 
the marginal cost of incremental power, an optimal curve may 
be drawn which is tangent to these "knees". | 

The example in figure 6-16 is for an energy storage of 
30 hours, that is, the storage system is capable of delivering 
power at rated capacity for 30 hours. Figure 6-17 displays a 


system of these curves for various levels of energy storage. 


6.7 Optimal level of energy storage 


The curves in figure 6-17 were reexpressed to plot 
demand as a function of energy storage with lines of constant 
availability in figure 6-18. The line of maximum value of 
energy storage revresents the highest output per unit of 
Storage possible for any availability/demand combination. The 
graph shows that for most levels of availability the maximum 
POwer output per unit cost of energy storage ranges from 3 to 


7 hours of storage. 
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Figure 6-16: Sample Storage Curve Construction. 
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Design Lanes For Storage Can- 
acity (Hours of Storage Required For 
Given Level of Demand). 
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Figure 6-18: Demand Fraction and Availability 


Tradeoff Lanes. 
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6.8 Variability of wave power (ore) 


The design lanes are based on the exponential power 
function, whose standard deviation (0,,) equals the mean power 
(P), or e/ P = 1.0. The tables cited in Chapter 4 on wave 
power availability indicate that ce SelnOner ores ts deep 
ocean and some coastal locations. However, there are some 
locations where the value of this parameter may be less than 
1.0. For this reason an adjustment factor needs to be applied 
to the design lanes to account for this change in power 
wariability. 

Using the simulation model the adjustment factor (K_) 
was determined and it is graphed in figure,6-19. The procedure 
for using the curve is as follows: 

1. Determine energy storage requirements from previous 

design lanes; 

2. Enter figure 6-19 with the actual value of Sa for 

the desired location and retrieve Koi and 

oe Multiply Ko by the energy storage from the design 


lanes to find the actual energy storage required. 
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Figure 6-19: Storage Reduction Factor (K_) vs 
ariability of Wave Power. 
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CHAPTER 7 


PEAS bo be ceRADEOEE STUDIES 


The objectives of any feasibility study include: 

1. To define a set of feasible alternatives which 

meet initial performance requirements; 

Tlomcaccermine sCOst-eftectiveness of the alternatives; 

and 

3. To assume definition of the alternatives to the 

level required for first order cost estimates. 

The feasibility design should be balanced. It is 
inappropriate at the. feasibility level to select specific 
subsystem components unless the user specifies such a selection 
as a deSign requirement. Subsystems should be identified by 
their general performance characteristics. 

The subsystems analyzed in this chapter include 
structure, conversion, generation, transmisSion and energy 
Storage. The effect of seasons On mean power available and 


a process for power rating selection are also discussed. 


yl Seasonality effect on selection of mean power level 


The design lanes in Chapter 6 are all based on a mean 


power level (P). This mean changes from season to season. 


Using figure 4-1 and assigning pa cic = 1.0, a typical 


seasonal breakdown may be described as follows: 





— 
— 





yes, 


season BCCAISOn 

Summer -60 Eee 
chien wane 1.40 Dead 
Fall/Spring 1.00 Pema? 


The question arises concerning selection of the design 
Pp, that is, should the system be designed for the minimum P 
(summer), the maximum P (winter), or the mean P (fall/spring)? 
Designing for the maximum P yields the most power delivered, 
but designing for the minimum P costs the least. ; 

The lowest $/KW system depends on the tradeoff between 
delivered power and soumenedloin//gneakzony/ aces (dem 
costs. The tradeoff is independent of structural costs 
because maximum power output diameter selection is independent 
of seasons. Energy storage costs are assumed to be constant 
Over the seasons. 

The mean power delivered varies with season. If the 
system is designed for the minimum P, then during summer months 
it is rated at 100% of the power available, during spring and 
fall it is rated at 60% of the power available, and during 
Winter it is rated at 43% of the power available. The compon- 
ents of this system will have a higher utilization than the 


components of a system designed for the maximum P (winter). 
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Besides this variation in mean deliverable power, the 
availability of the power delivered changes with season for a 
given design. For example, if the system is designed for the 
minimum P and a value of X is assigned to the availability of 
the rated power during summer, then during winter months the 
availability of the rated power will be greater than X because 
the mean power level in the ocean is higher. 

Assuming that the marginal value of availability 
equals the marginal cost of power, then the equivalent power 
delivered for each of three designs (min P, mean P, and max P) 
may be adjusted for availability by using figure 6-18 and 
holding energy storage constant. 

Table 7-1 lists the power rating, the deliverable 
power, the availability adjustment factor and the equivalent 


power delivered for each of the three designs (min P, mean P, 


and max P) over all seasons, non-dimensionalized to the mean 


P design. 
TABLE 7-1 
RELATIVE POWER DELIVERED FROM SEASONAL DESIGNS 
Equivalent 
Delivered Availability Delivered 

Design Power Rating Power Comme ce1On Power 
Max P 1.40 Te ditt no ioe 
Mean P 1.00 1.00 1.00 in O0 


Min P .60 .67 1.25 .83 
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If conversion/generation/transmission (CGT) costs 
represent 0% of the total system cost then the max P design 
is optimal since it delivers the highest equivalent power. 

If CGT costs are 100% of the total system cost then the min P 
design 1S optimal since it has the least cost. 

Using data in table 7-1, the costs of power for CGT 
costs representing varying percentages of total cost were 
calculated and displayed in table 7-2. The minimum cost/KW 
designs are circled in the table. Since CGT costs are expected 
to lie between 10% and 40% of total system cost the selection 
of the annual mean power level for design considerations is 


the optimal selection. 


TABLE 7/7~2 


RELATIVE $/KW FOR COMPETING DESIGNS 
CGT costs as % of total cost 


Design O 10 20 30 40 50 
Max P 1.02 1.06 1.09 eS lpi: 
Mean P 1.00 1.00 


Min P e201. 1 6 ie] On? 1.02 


oe, 
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7.2 Seructure 

The most critical structural design element fOlee ene 
WPG will be longitudinal bending moments. The Load Line 
Regulations (ref. 4) assert that good design practice requires 
that values of length/depth lie between 10 - 13.5 for ocean- 
Meeng vessels. For vessels operating in the Great Lakes 
regions where wave heights are limited, the Regulations 
require L/D less than 13.5 for vessels less than 325 feet in 
length and L/D less than 19 for vessels with lengths greater 
than 600 feet. One of the reasons for these requirements is 
that insufficient depth lacks the necessary moment of inertia 
to give proper Lometbewdltnael stiffness. 

The strength of the cam lies in its circular cylindrical 
backbone. Ocean-going vessels midship cross-sections are more 
nearly approximated by box shapes than circles, implying that 
for similarly sized structures the cam will have a lower mid- 
ship section moment of inertia than the ocean-going vessel. 
With this in mind it would therefore be prudent to restrict 
the cam's maximum value of L/D to 10. If interconnections 
between cams are deemed feasible then the total L/D rather 
than the individual cam's L/D must be less than 10. If cam 
interconnections are flexible then the L/D requirement may be 


relaxed since the flexible link will not carry a load. 











weight 


Is 


Structural costs are frequently based on structural 


(W). A set of cost relationships based on Maritime 


Administration data (ref. 17) for ocean-going tanker design 


in 1977 follows: 


Material Cost M1 
Direct Labor Cost M4 
Indirect Material Cost M5 
Indirect Labor Cost M6 
Material Engineering Cost M7 


Overhead Cost M9 


= 530W 

= 384W 

= 31.8W 
= 69.12W 


= 10.6W 


= 45.31W 


motal Cost sue 


= 1070.83W 


Structural weight (W) is estimated from the following 


Pande Olen t/ Dyrlal2—. O63. £7 D) 


Coorg )) (4D) 


equation: 
w= vcr’°> 7°© 4 1008) 
where Y = 3.8 - 1.1(L/D) + .1(L/D)* 
aa eo yee 
F = 3.8 + 2.1 L/1000 
2 = mos L*D(C, +.7) ime 


(L is in feet) 
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For the cam design, Cy =—eewaoe folding L/D constant 


at 10, the equation reduces to 


W= 1.075 x 107° (14¢2:24) (3.8 + feat) °° (12D) °° 
i 


L 





(L) 
(L & Din feet) 


This equation was computed for a range of diameters 
from 1 ~ 50 meters with the resulting weight/length a function 
Bit D 


Weaght (tons) _ 2 
Length (m) = .1865 D Peper mecers) 


The proportionality constant, .1865, represents the 
Bverage of the constants calculated for the various diameters. 
The standard deviation of the constants calculated was +3% of 
the mean, verifying the proportionality of structural weight 

2 


myth D for the cam. 


The total cost/length for the cam structure is, 


ALS $1070.83 W 


on Ome (er hei6 5 D*) 


Sl ke ee ik D* per meter of length 
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7.3 Conversion 

The characteristics of the conversion system are 
governed by the following design Peon anenkes: 

(1) the cam's rocking motion must be converted to 
rotational motion at speeds compatible with 
generator requirements; 

(2) damping characterized by variable damping coeffi- 
cients must be utilized to achieve high 
efficiencies; 

(3) the system must operate with variable speed inputs 
characterized by @a/H = 1.0 for maximum efficiency; 

(4) the system must control power surges; and 

(5) the system must operate in an ocean environment. 

Only one conversion method, hydraulic/pneumatic 
Satisfies all these requirements. 

In order to fulfill requirement (1), variable speed 
radial piston, axial piston or vane pumps may be used. The 
Slow speed input implies ee piston pumps will be best 
Suited to the conversion. | 

The pump power takeoff will require a fixed reference 
which may be fulfilled by floating an object behind the cam 
or submerging a ballast weight beneath the cam. In either 
case the fixed reference will be attached to the ends of the 
cam and through seals and bearings the cam will rotate about 


the reference. 


‘ 
_ 
yy 7 
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Hydrostatic drives are best suited to applications 
requiring wide ranges of speeds and torques. The basic 
elements of the hydrostatic drive are a pump, : motor and a 
pump driver. Mechanical energy from the cam driver is con- 
verted to fluid energy by the pump and reconverted to 
mechanical energy by the motor for use by the electrical 
generator. Since it is possible to provide infinitely 
variable pump and motor displacements, speed ratios may be 
stepless. 

There are four basic pump/motor configurations: 

(constant capacity pump, constant capacity motor: 

This configuration provides a fixed ratio drive speed 
which is proportional to the displacements of the pump and 
motor. If the two units have equal displacement then the 
drive ratio is one and output speed equals input speed. Speed 
variation may be accomplished by using bypass flow in which 
case power is diverted away from the motor and wasted in the 
form of heat. 

(2) Variable capacity pump, constant capacity motor: 

This configuration provides clutching, constant torque 
and variable speed. With the pump at zero displacement, the 
System is clutched and motor speed is zero. An increase in 
pump displacement increases motor speed, up to the design 


limit, with approximately constant torque. 
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(3) Constant capacity pump, variable capacity motor: 

This system delivers constant power over the design 
speed range. Reduction in motor displacement results in 
increased speed but reduced torque. 

(4) Variable capacity pump, variable capacity motor: 

This arrangement is particularly well-suited for 
designs requiring a wide range of speeds and torques. It 
combines the constant power characteristics of system (3) and 
the constant torque characteristics of system (2). 

Of these configurations system (4) is best suited to 
the WPG. The variable pump allows selection of the optimum 
damping ratio and the variable motor allows operation over a 
wide range of input speeds without wasting power. Use of this 
system implies coupling with a variable speed generator. 

For the constant speed generator application system (2) 
appears best suited. The variable pump could be controlled 
to select optimum damping. A valving mechanism for wasting 
excess energy, or an energy storage link, would be necessary 
because of the constant motor speed requirement. 

In any configuration gearing is possible between the 
fluid motor and the electrical generator. 

Table 7-3 lists the characteristics of hydrostatic 


conversion. 
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Direct mechanical conversion will not be feasible 
because the cam's velocity reaches zero twice with each 
passing wave. Ina purely mechanical system the rocking 
motion could be converted to rotational motion with an 
arrangement of ratchet/pawl and worm/crown gear mechanisms. 
Mechanical, infinitely variable speed transmissions, like the 
Voight unit utilized in wind power, are capable of applying 
damping when coupled with field modulated electrical generators. 
But with zero velocity twice each cycle a clutch would be 
needed between the drive and the generator to account for 
inertia effects. Since the transmission must be non-slip in 
order to apply appropriate damping, the clutch requirement 


causes the purely mechanical conversion system to be infeasible. 


7.4 Electrical Generation 
Acceptable schemes for generation compatible with power 
grids may be generally classified as constant speed constant 


frequency (CSCF) and variable speed constant frequency (VSCF). 


7.4.1 CSCF Systems. Synchronous and induction 
generators belong to this category. The constant speed 
requirement dictates the use of energy storage between the 
cam's rocking motion and the generator shaft or a valving 
mechanism which wastes excess energy. Suitable control systems 


are required for sensing and reacting to input speed variations. 
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7.4.2 VSCF Systems. There are several types of 
Variable speed machines in existence today, all of which are 
more costly than the CSCF machines. 

(1) AC-DC-AC. The AC-DC-AC conversion system 1s an 
AC generator with DC field excitation. The variable frequency 
DC output is rectified and converted to AC through the use of 
an inverter or a motor/generator set. The system rapidly 
reaches stability with the variable speed input. AC-DC-AC 
may be attractive in deep-ocean installations where DC trans- 
mission to shore is most economical. 

(2) Double Output Induction Generator. This machine 
makes use of an induction generator with rotor energy regen- 
eration. Both rotor and stator supply power to the electrical 
grid. Power from the stator which is generated at grid fre- 
quency is directly fed to the grid. Power from the rotor 
which is generated at slip frequencies is first converted to 
DC, then inverted to AC before being fed to the grid. 
Efficiencies greater than 80% can be achieved with input speed 
Variations from 20% to 100% of rated speed. 

(3) Pole-changing Generator. This generator electro- 
Magnetically alters the field pole pattern according to shaft 
speed. The electrical output is constant in frequency. There 


are no restrictions on input speeds. 
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(4) Field-modulated Generating Systems. These units 
are generally excited by an AC power source fixing the output 
frequency and modulated by rotational frequency. High shaft 
speeds are required. Output efficiencies are comparable 1o 
synchronous machines and constant over input speed ratios 
me o:1. 

Comparisons of off-speed efficiencies for several 
generators are displayed in figure 7-1. 

Table 7-4 compares costs and efficiencies for a number 


of generation units. 


7.5 Optimal Power Rating 


It was shown in section 7.1 that the WPG will deliver 
the most power over all seasons if the design is based on the 
mean annual power level. But the selected design power level 
may be more narrowly defined than this. 

Power rating can be optimized. If the power rating 
selected is a small fraction of the mean annual power available 
(P) then the WPG will deliver very little power. If the power 
rating is selected as several times P then the system will 
deliver more power but at a higher cost. For the analysis 
that follows the $/KW of the conversion/generation/transmission 


(CGT) subsystem is assumed to be constant. 





Efficiency 
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Figure 7-1: Off-Design Performance of Several 


Blectrical Generators. 
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The fraction of total power delivered up to any given 
level of P may be derived from the simulated Zero storage line 
in figure 6-17. Expressing power rating as PR = fraction of 
power available, the resulting curve is plotted as the solid 
line in figure 7-2. 

The two dashed lines in this figure are included for 
comparison of simulated results with the theoretically derived 
exponential power distribution. The upper dashed line repre- 
sents the fraction of all wave arrivals for any given PR from 
the equation, Fraction = l - caus (see Chapter 3). The lower 
Seened line represents the theoretical fraction of power in | 
those waves for any PR. The solid line which also represents 
the fraction of power in waves lies above the lower dashed 
line because of the filtering characteristics of the WPG 
efficiency curve. An empirical formula for the solid line 
was derived: 


eR 


ee (l= 2 a 


where Y = fraction of total power available 


K = .90 





Fraction of Power Available 
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Simulated Fraction Of Power 
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Figure 7-2: Fraction of Power Available in Waves uv 
to a Certain Level of Power Rating (PR). 
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Once the power rating is selected the fraction of total 
power captured may be calculated. The following assumptions 
make the calculations more tractable: 

(1) the CGT subsystem delivers all the power in waves 

having power levels less than rated capacity; 

(2) for waves having power levels above the rated 
capacity, the system delivers the rated capacity; 
and 

(3) there is a maximum power level above which the cam 
ceases to operate. 

Assumption (3) requires some clarification. The WPG 
has not yet been tested for operation in severe sea states, 
but it is reasonable to expect that there is some maximum sea 
condition above which the cam cannot operate. In this situa- 
tion the cam may be submerged or allowed to float with power 
generation apparatus secured. This maximum condition may be 
expressed as PRax’ The fraction of power above ER is. lest. 
Let Y = fraction of total power available which is 

captured by the CGT subsystem 

PR = design power rating 

PRiax = maximum power rating 

¥Y = (1) all of the power up to rated capacity + 

+ (2) some fraction of the power between rated 


capacity and maximum capacity + 





gl 


+ (3) none of the power above the maximum capacity. 


Term (1) above = K(1 - saa 


Term (3) = 0 
Term (2) was derived in its general form. Referring to 
figure 7-3, term (2) graphically is the ratio of Area A/Area 


(A+B). In numerical form, term (2) is 


R PR 2 


K(PR) fe * - e@° max] 


-PR -PR 
(l -e max) (1 + PRuax> + PR-PR ax” (l-e ) (1+PR) 


Therefore, 


eek PR 2 
a Ge PR) -e max] 


: (l-e max) (1+) "max)+PR-PR__|-(1-e ~ *) (1+PR) 


(1) 
Although the value of PR is not known, PR = 3 1s 
max max 
a reasonable assumption. Referring to figure 7-2, eRe es 3 
corresponds to 95% of all waves and 84% of all power. [In 
Other words the assumption of Reet 3 implies that the 
System will be completely shut down for 5% of the time due to 


excesSive sea conditions and 16% of the total power available 


Will be lost. 
Equation 1 was calculated for ee = 3 and plotted in 
figure 7-4. The dashed line in this figure represents 
PR 


Y = Pe en ), which corresponds to a system which is shut 


down when power levels exceed PR. 
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Since COsStcs are assumea to be Linear with size for the 


CGT subsystem, the slope of the "Y" line in figure 7-4, 


adYy/dPR, is 


actually the marginal cost of the CGT subsystem, 


_ $/KW cgt 


MP ogt = S§7KW system - cgt oe 


c-— = Ke + 


where A= 


= K(e -e max 


-PR Ves) [Pe el eu 


(2) 
aie 


PR 


eee ee oR rrprePri = (lees) (1l4+PR) 
max max 


a 


JEN 5 SRNR SENS 


(e -e max) 


Teo aS ) *_2KPRe- 


gue be Z 


= RE e -e max) 


-pre *R 


Values of dY/dPR (eq.2) have been plotted in figure 7-5. 


Example use of figure 7-5: 


as) 
(2) 


(3) 


Calculate MF : 
(ele 6 cs 


With MP oot enter the vertical axis and mark the 
corresponding point on the MOE ae line. 
Read down to horizontal axis for the design power 


rating, which is just PR x P. 
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Figure 7-5: Design Lane For Selecting Optimal 
Power Rating And Bandwidth Efficiency Given 
Marginal Cost of CGT Subsysten. 
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(4) Mark the point on the "Y" line which is vertically 
above or below the point on the BY ets line and 
read across to the value of Y. 


Numerical example: 


MP ogt = .16 


Design power rating = 1.18 P 


Fraction of power collected (Y) = .77 


The design power rating is optimal for the assumptions 
stated previously. This number multiplied by the cost/KW of 
the CGT subsystem yields the actual cost of the subsystem. 

The fraction of power collected, Y, may be regarded as 
a subsystem efficiency, and it will be referred to in Chapter 
8 as the bandwidth efficiency, Tears 

Efficiencies of specific generation, conversion and 
transmission equipments are assumed to be 1.0 in this analysis 
and they must, therefore, be taken into account in the final 
design process. 

The preceding analysis may be utilized to show that 
variable speed generators should be used in the WPG system. 
A constant speed generator will operate between a cut-in and 
a cut-out wave power level. Below the cut-in level no power 
is delivered. Between the cut-in and cut-out levels the 
‘generator produces constant power equaling the cut-in power 


level. 
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Numerically, the amount of delivered power for the 
constant speed generator without energy storage is 


PR 


¥-K(l-e  )= 


K(PR) (e-P 8-e max) 7 


PR 


= gle 
(l-e max) (1+PR_ |) +PR-PR, x (ine yee R) 


a 


Graphically, the constant speed generator's delivered 
power corresponds to the difference between the "Y" line and 
the Matt-e* *)« line in figure 7-4. Figure 7-6 compares the 
plots of the variable speed generator (Y) with the constant 
speed generator (Y - KGieee ) Clearly, the constant speed 
generator without energy storage cannot compete with the 
variable Bee nera tee. In addition, the spilled eens 
with the constant speed generator must be dissipated as heat 
in the hydraulic conversion system which increases conversion 
costs. 

It was pointed out in section 7.3 that energy storage 
for the constant speed generator must include a hydraulic link. 
But hydraulic energy storage is prohibitive in cost (see table 
7-6). If the hydraulic conversion link is coupled with 
kinetic energy storage, like a flywheel, then the flywheel 
Must physically be located in the cam. The input rotational 
speeds of about 1 RPM from the cam to the hydraulic link would 


have to be increased to about 20,000 RPM for efficient flywheel 


Operation. This is infeasible. 


Fraction of Power Delivered 
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Variable Sveed Generator 


——--—— Constant Sveed Generator 
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re 7-6: Comparison of Power Delivered by Var- 
iable Speed and Constant Speed Generators. 
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With variable speed generators the flywheel, or any 
other energy storage device, could be located ashore and many 
WPG units could feed a single storage device. Flywheel speeds 
could be controlled electrically rather than mechanically. 

For these reasons the variable speed generator is 
preferrable to the constant speed generator. From table 7-4 
the brushless DC generator with processing for inversion to 


grid AC appears to be the most feasible generation subsystem. 


y.6 Transmission 

The transmission subsystem depends on system configura- 
tion. For distances less than 20 miles, AC electrical trans- 
mission is the most cost effective method. Beyond 20 miles 
from shore, DC transmission is most cost effective. 

It is not necessary to require transmission of 
electricity. Power generated at the cam may be used in 
hydrolysis of water for production of gaseous hydrogen (GH.) 
or liquid hydrogen (LH. ) one sues some derivative of hydrogen 
like ammonia. GH. may be feasibly piped over long distances 
and LH. may be transported by barge. The hydrogen produced 
may be sold commercially or used to run a fuel cell for 
electrical generation. 


Characteristics of the various transmission schemes 


are displayed in table 7-5. 





—_ 
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y./ Energy Storage 

Since the power in ocean waves varies SO considerably 
with time, some form of energy storage may be required in the 
system. As previously mentioned the CSCF design with hydraulic 
pressure energy storage is infeasible. In VSCF systems energy 
storage may be accomplished by means of flywheels, batteries 
or containerized hydrogen (GH. Or LH). 

Comparisons of various energy storage methods are 
outlined in table 7-6. 

For GH. storage additional costs are incurred in 
electrolysis and reconversion to electricity. For LH. storage 
additional costs are incurred in electrolysis, liquefaction 
and reconversion. These added costs are outlined 4 wal Die wy oi 

From table 7-6 the three most likely forms of energy 
storage are lead-acid batteries, modern flywheels and LH.. 
LH. is included here because it may become the optimal storage 
medium if the WPG is sited at great distances from shore. 

The technology of lead-acid batteries is well known. 
The only major limitation imposed by the use of these 
batteries is the slow charging rates (in hours for the costs 
associated with batteries in table 7-6). lLead-acid batteries 


with charging rates on the order of minutes do exist, but the 


costs are several times the cost listed in table 7-6. Slow 
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sEfea 


charge batteries may be feasible if electrical inputs come 
from a large number of sources scattered over a large geo- 
graphic area, which may tend to smooth out the power delivery. 

A considerable amount of interest has been generated 
in the past few years over the use of high energy density 
flywheels. Energy densities have been dramatically increased 
through the use of composites like fused silica and protection 
costs have been reduced because at failure stresses the com- 
posites shatter in millions of small pieces instead of large 
steel chunks like flywheels of the past. Run down times for 
composite flywheels in vacuum chambers have been shown to be 
on the order of months. The limitations on flywheel energy 
meolage are: 

(1) very high rotational speeds (20,000-25,000 RPM); 

(2) vacuum chamber sealing; and 

(3) inherent rotational inertia. 

The high rotational speeds can be achieved electrically 
by using a multipole generator to drive a permanent magnet 
motor in the flywheel. The high generator/motor pole ratio 
achieves the desired high rotational speeds. The interface 
is electrical, reducing the vacuum chamber sealing complexity. 
During charge the flywheel acts as a motor, and during dis- 


charge it acts as a generator. The use of the multipole 
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generator is an added cost to cost of flywheel storage. The 
rotational inertia limitation can be eliminated by using 
counter-rotating flywheels or by rigidly securing the unit 
to the earth. 

Taking into consideration all capabilities, limitations, 
costs and efficiencies for the various energy storage schemes, 
the modern flywheel appears to be the most viable system for 


the WPG. 


7.8 Other subsystems 

Mooring and direction-keeping subsystems cannot be 
analyzed until free-floating tests are conducted on the cam. 
In Ocean Thermal Energy Conversion (OTEC) systems, mooring 
costs typically run about 10% of structural costs; this cost 
will be assumed for the WPG. 

There will undoubtedly be a reduction in the effective 
mean power density at a particular site due to the inability 
of the cam to maintain itself parallel to wave crests and its 
inability to maintain a perfect fixed reference, assuming the 
fixed reference is not a structure fixed to the ocean floor. 
These effects, however, will not appreciably distort the 
Overall design. If the power density is effectively reduced 
then all system components, except structure, will be down- 


sized accordingly. 
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TAbias 7 =3 


CONVERSION SYSTEMS (no storage) 


S/KW Efficiency 
Pneumatic (es 5 
Oil hydraulic 64 ASS; 
TABLE 7-4 


ELECTRICAL GENERATION SYSTEMS 


$/KW ELeacs ene, 
Synchronous 60 ~94 
Induction 39 28s 
FMGS 72 80 
Pole changing 900 150, 
DC/AC P73 . 86 
Multipole 86 {OS 


Rectifier 30 ~92 
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TABLE. 7 =5 


TRANSMISSION (100 MW nominal sizing) 


DC electric cable 
GH. pipeline 


LH. barge 


2 
Fuel cells 
Electrolyzer 


Liquefaction 


Factor increase in cost for designs 


1 MW: 
50 MW: 


500 MW: 


4.3 
1.4 


sii 


S/KW 1 mile 


nig 


30 


IO 1 
oal 
114 
9 

89 


234 


$/KW 1090 miles Ey ee ene 


a 
sks 
~95 
70 
84 


80 


other than 100 MW: 





\ 





Lead-acid battery 
Modern flywheel 

LH. 
Sal hydraulic 


Compressed air 


GH, 
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TABLES 7-6 


ENERGY STORAGE 


$/KW = Lb/KWh 
37 67 
80 20 
82 06 
490 dg 
674 143 
7 35 2.9 


Ft?/KWh 


tee 


‘asm 


67 


eye 


08 


Be lveciency 


60 
225 
oD 
515300, 
20 


40 
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CHAPTER 8 


COST ANALYSIS 


Bre L Costs 

Capital costs include costs for structure, conversion, 
generation, transmission and storage. Annual costs may be 
expressed as a % of capital costs. Data was assimilated from 
various sources for ocean~based power generating systems to 


meeodauce the following annual costs: 


some caplcaleoste 


Taxes 2D 
Insurance fe 18) 
Operating and Maintenance > .5 
Depreciation and Overhead 20 
Return 11.0 

237.0 


The annual fixed charge rate of 23% must be adjusted, 
or annuitized, for the life cycle of the investment. The 


annuitized fixed charge rate (AFCR) = 








gO! 


where 1 fixed charge rate 


3. 
II 


investment life (years) 


The investment life for the WPG is assumed to be 


20 years which yields AFCR = .2337. 


B.2 Power Analysis 

The WPG may be analyzed by two methods: (1) fuel 
saver credit - a conventional plant is shut down with the 
resulting savings in fuel; and (2) baseload capacity credit - 


a-conventional plant is physically replaced by the WPG. 


8.2.1 Fuel saver credit analysis. The simplest 
application of the WPG is fuel saving. Power from the WPG 


is fed directly into the utility grid and depending on the 
grid load, conventional fuel-fired plants are backed off with 
a savings in fuel resulting. The value of the fuel savings 
to the utility network requires knowledge of the Incremental 
Heat Rate (IHR: BTU/KWh). The IHR represents the number of 
BTU's of fuel saved with the reduction of one KWh of conven- 
tional generation. IHR varies constantly. It is currently 
about 9500 BTU/KWh for off-peak power, 11000 BTU/KWh mean 


power and 12300 BTU/KWh for peak power. 
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The value of any fuel-less power generator as a fuel 


saver 1s 


Fuel Saver Value (¢/KWh) = [(Fuel Cost (¢/10° Ban inn 


THR (BTU/KWh) (1) 
6 
10 


Figure 8-l shows this relationship between fuel costs 
and fuel saver values for various IHR's. 
The busbar cost equation is the basic formula for all 


power cost analyses: 


_ {Capital Cost ($/KW) xAFCR(Fraction/yr) ]xl00¢/$ 
poe cost (*/KWh)- availability (Fraction)] = 8760 hrs/vr 


(2) 
To simplify the use of the equation, the following 


abbreviations will be used: 


FSV = Fuel saver value (¢/KWh) 
BLEC = Baseload energy cost (¢/KWh) 


AV 


Availability (fraction) 


ce 


Capital cost ($/KW) 


AFCR = Annuitized fixed charge rate (fraction/yr) 
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Figure 9-1: 


Cost of Oil ($/SBL) 


System Incremental Heat Rate 
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In fuel saver analysis the following variation of 


the busbar equation applies: 


CC (AFCR) 


FSV = TAO ob G3) 


Once the fuel saver value (FSV) is known, the 
corresponding breakeven fuel cost may be determined from 


figure 8-l. 


8.2.2 Baseload capacity credit. If the availability 
of power from the WPG is sufficiently high then base load 


credit may be applied. The pertinent .equation is, 


_ CC (AFCR) 
BLEC = 3776 (AV) (4) 
8.3 Availability 


At first glance there is no difference between 
equations (3) and (4). The actual difference lies in the 
rather complex notion of a€vailability as it applies to the 

mG. In fuel saver analysis it does not matter when the power 
from the WPG is introduced to the grid: as power arrives from 
the WPG, conventional plants are shut down. For this reason 
the availability of the mean power level in the ocean is 
equal to 1.0 for the fuel saver analysis (neglecting adjustment 


factors). 








eet 


Therefore, AV es = Ae. 

In the baseload analysis the concept of availability 
for the WPG changes: it is the level of power which the WPG 
can be relied upon to deliver for a certain fraction of the 
year. The level of demand achievable and the energy storage 
required for a given value of availability may be extracted 
from figure 6-18. 

It was mentioned in section 7.1 that the mean power 
level changes with seasons. For this reason baseload credit 
may Only be applied to the lowest seasonal value, that is, 
to .60 P, -the level in summer months. Remaining power may 
be assigned fuel saver credit. The amount applied to baseload 
is .60 (A) (DF), where A = availability from figure 6-18 and 
DF = demand fraction from figure 6-18. The amount applied to 
fuel savings is (1-AFCR) (1-.60(A) (DF)), where (1-AFCR) adjusts 
for the difference between fuel savings and baseload value. 


The availability in the baseload analysis then is, AV 


bl 
po0(A) (DF)+(1-AFCR) (1-.60(A) (DF)). 


There are two adjustment factors which must be applied 
to both definitions of availability: 

(1) A certain fraction of the year must be set aside 
for equipment maintenance and down-time. Conventional plants 
which are mechanically more complex than the WPG assign 


AV = .75. Ocean Thermal Energy Conversion plants which are 
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less complex mechanically than the WPG assume AV = .95. It 
seems reasonable to set aside a mechanical availability Om 
moo tor the WPG. 

(2) In Chapter 7 a method was introduced for optimizing 
power rating (figure 7-5). The bandwidth efficiency, ana" is 
an adjustment factor which must be applied to availability. 

The two adjustment factors may now be applied to AV 6. 


and AV,, (the bars over Ae anda AVeeainethe following 


b fs bil 


equations imply adjusted versions). 


AV 5. = .85 "hand (5) 


ou () (DEy erat eR) (l—.60 (A) (DF) ) 


(6) 


AV = .85 eana 
The value of (A) in equation (6) for use in figure 
6-18 is .75, since baseload analysis of the WPG should be 
conducted with respect to conventional plants which assume 
oy = .75. 
Substituting equations (5) and (6) into equations (3) 


and (4) yields the relationships for WPG cost analysis: 


Sakon AER) (7) 
87.6 (AV,.) 
pec = 6C_(AFCR) 
87.6 (AV,,) 
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8.4 System Design 


8.4.1 Site selection. The first step in the system 
design process is to select the site for operation of the WPG 
from the tables in Chapter 4. Data should be extracted from 
the "All Seasons" table since the WPG delivers the most power 


when designed to mean annual power levels. 


8.4.2 Diameter selection - without separation effect. 


occurs at wva/g = .64. 


Utilizing the dispersion relation, D a = .,2034 Tr. The 


In Chapter 6 it was shown that Dias 
minimum cost diameter was shown to be 60% of Dass delivering 
68% of the power of De ass (figure 6-5). Therefore, design 
Memeter me .122 7 and the power available = 68% of the 
maximum power. This delivered power may be termed a system 
design efficiency and will be defined lavectem: Therefore, 
Ndesign ~ OEE 

Structural cost = $199.71 oan Applying the factor 
increase for mooring cost, structural cost = 1.1(199.71)*x 
(6122 T°)? = 3.27 t’/p 


where P = power delivered to 


del’ del 


user. 
The cam efficiency is extracted from figure 6-8 using 


the site selection parameter O/T. 


— 


- +> 
—_— > 


> 





Pn 
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8.4.3 Diameter selection - with separation effect. 


With separation, Dee (3-D) = .68 2 eee (2Z—)))7 and the minimum 

cost diameter occurs at .52 De ee delivering 56% of the power 

at D - Therefore D (3-D) = .68 (.2034 7) =e og 3 Te, 
max max 


and the design diameter = .52 (.1383 T7) = .072 T*. 


Structural cost with separation = 1.1 (199.71).« 
2.2 4 

“@c072 Ta)= = 1.139 Ff /P so] 

The cam's efficiency is extracted from figure 6-1l 


using the site selection parameter O/T. 


8.4.4 Conversion. The hydrostatic drive is the 
selected conversion system with a cost of 64 S$/KW and 


efficiency of .65. 


8.4.5 Generation. The DC/AC generation scheme is the 
selected method of electrical generation with a cost of 133 


S/KW and an efficiency of .86. 


G.450 POwer rating selection. Power rating (PR) is 


obtained from figure 7-5. The bandwidth of useable power 


Mana? is also determined from figure 7-5. 
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8.4.7 Transmission. DC Cable transmission is the 
selected method having a cost of 101 S/KW with an 


mritrerency of .95. 


8.4.8 Storage. The flywheel storage system is the 
most feasible choice for the WPG. Storage is not necessary 
for the fuel saver analysis, but it may be necessary (or 
desireable) in the baseload analysis. The cost of storage is 
80 S$/KWh with an efficiency of .95. The actual storage quantity 
in hours is determined by multiplying the number of hours deter- 
mined from figure 6-18 by the storage reduction factor (Ko) aL)9) 
figure 6-19 using the site selection parameter O,/P. The fly- 
wheel input/output device is the multipole m/g set having a 


cost of 86 S/KW with an efficiency of .95. 


B.5 Design Method for Fuel Saver Analysis 
8.5.1 With separation effect. 


"design Oe 
“conv = ge 
gen = .86 
Ntrans ~ °?>? 
i : from figure 6-11 


Noverall “design "cam "conv "gen "trans 


a0 Nee 


Padel ~ "overall 
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: _ 1.139 7° 
structure eal 
S S ° 
trans gen conv 
= ———— nS 
ican ~ ‘Woe n n n ae 
rans trans gen erans, Gen Conv 
-car = 389(PR) 


Use figure 7-5 to determine PR and Dee a" 


ME oer _ —— : 
structure 
perform iteration with values of MF oer and PR to make this 
equation valid. 
ieee ee 
CC(S/KW) = a + 389(PR) 
del 
AV es 7 ols “band 
psv = SC(AFCR) 
87.6 (AV...) 
8.5.2 Without separation effect. The following 


changes are made to the analysis of section 8.5.1: 


design ee 
Noverall ~ 2a oam 
n from figure 6-8 


cam 





Jey 


4 
cc = 2:27 T 4 389 (pR) 


Pdel 


8.5.3 Example fuel saver calculation. North Atlantic, 


all seasons, with separation. 


8.5 sec 


KH 
ll 


td 
il 


37.1 KW/M 


O/T = .31 “ 


From figure 6-ll, for O/T =o 1s, eee oO 


ee, NoP = 9,24 KW/M 
: eee aes conn 
structure P 
del 
_ 389 
MFeer ~ 643 ‘PR? 


Performing iteration in figure 7-5 yields PR = .68, Te a 64. 





IES 


CC = S$ iryceure * Scop = 543 + 389(PR) = 908 $/KW 
AVe, = -85 nyoig = (-85)(.64) = .54 
Bey = CoURECR) — eee = 4.49 ¢/Kwh 

87.6 (AV, .) ees 


FSV = 4.49 ¢/KWh 


The corresponding breakeven fuel cost for fuel savings 
= 4.49 ¢/KWh may be determined from figure 8-1. For IHR = 
11000 BTU/KWh the breakeven fuel cost is 4.05 $/10° BTU, or 


Mono, sei Of oil. 


8.6 Design Method for Baseload Analysis 


This method is the same as that for the fuel saver 
except that AV) must be calculated uSing equation 6. Costs 
must be calculated for zero storage and for the highest 


return storage (figure 6-18). The lower cost design should 


be selected. 


8.6.1 Example baseload calculation. North Atlantic, 


all seasons, with separation: 





Eo 


(l) No storage 


T = 8.5 sec 

P = 37.1 KW/M 
O/T = .3l 

o,/P = .98 

No = .25 

Pael = 9,24 KW/M 
Thana = .64 


CC = 908 $/KW 
AVE = 


from figure 6-18, for A = .75 and zero storage, 


85 lea Cusco mba y(t ARCR) (1-.60(A) (DF) ) 


DE ao. 


AVL 43 


CC (AFCR) 
AV, , ) 


87.6(A 


BLEC = 
lop 


5.63 ¢/KWh 


(2) With storage 
From figure 6-18 the storage quantity yielding 
Maximum return is 5.26 hours, at A = .75 and DF = .67. From 
figure 6-19 with o,/P = .98, K. = .982. Actual storage 


S 
required = K (5.26) = 5.16 hours. 
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Storage cost = 80 S$/KWh 


"storage y 242 


Storage I/0 cost = 86 S$/KW 
Nt/0 = .95 


60% of power is diverted to storage: 


— (.6) (86) = 
eens | ete (2) (86)! = 95 $/KW 


oc NeoS) 


AV, 4 = .85 nong [-60(A) (DF) +(1-AFCR) (1-.60(A) (DF) ]=.46 


2 
Ole -6(.30) “eam “storage "t/0 t .4(.30) “eam 
= 6 A) 
Pea igee eee) KW/M 
s25 _ 
cc = 908 (55) + 426 + 95 = 1514 S/KW 
peg) -o<BECR) = 8.78 ¢/KWh 
87.6(A bL? 


The cost of baseload energy is lower for the system 
without energy storage than for the system with energy storage 
implying that even the best energy storage scheme, the fly- 

wheel, 1s still too costly for economical use. It is worth 


| 


| 
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noting, however, that the demand fraction (DF) without energy 
storage is .3l, while DF with energy storage is .67. fThis 
means that in order to achieve a certain level of baseload 
credit, more than twice as many no-storage WPG's would have 
to be assigned to the baseload as would be required with 


storage-configured WPG's. 


ee f EPeegenetical Ocean Energy Warm Cost Calculation 


The energy farm electrolyzes water to produce LH, OG 
shipment to shore and conversion to electricity in fuel cells. 
The analysis is conducted for fuel saver with separation 


effect and without energy storage. 


Assumed ocean parameters: 


T = 8.5 sec 

P = 100 KW/M 

O,./T = 926 

opr Oc 

o,/P = kat 
ectrolveis = O° 7/7 kNa on = .84 
"liquefaction = 234 $/KW, n = .80 
eed Pedi 7 $/KW, n = .70 
° = 114 S$/KW, n = .95 


transmission 





| 





- 








PZ 


Ny = +13 
Pac, = 13 KW/M 
: ee 
oie, abvon eh bbate 13 
Sogn = 1668 (PR) 
_ 1668 | _ 
MF oer = TA (PR), yields PR = .34 


"band = 44 


CC = 714 + 1668 (.34) = 1281 $/KW 


AV e. =4s0 5 eee oi?) 


FSV = 9.14 ¢/KWh 


With this high cost occurring in such favorable 
ocean conditions it is unlikely that production of LH, will 
be feasible. Costs of electrolysis, liquefaction and fuel 


cells must be drastically reduced and efficiencies improved 


before LH. production will be feasible. 





= 
2 


—_ 
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B.8 Value of Energy Storage 


Using the calculations in the previous baseload 
credit analysis, the breakeven cost of energy storage may 
be computed. 


ce. + > storage : AVE (with storage) 
ce 


AV, 5 (without storage) 


= — > (908) - 908 = 63 $/KW 


~ 46 
~43 


> storage 


This is the breakeven cost for 5.16 hours of storage 
plus. the cost of the I/0 device. Assuming the I/0 device 
represents 20% of total storage cost, the breakeven value of 


One unit of storage is, 


63(.8) 


5.16. Sits 


It is unlikely that any energy storage device will 
achieve this low cost of 10 S$/KWh. For this reason energy 


storage is considered infeasible for the WPG. 


p8.9 Structural Cost Reductions 
The optimization methods for cam diameter in Chapter 
6 were based on the cam being built entirely of steel. Calcu- 


Mations of costs in the fuel saver mode showed that structural 
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costs are about 70% of total cost even with the separation 
effect. Structural costs must be reduced. 

It should be possible to build a cam having a small 
cylindrical steel strength member which houses machinery and 
a shroud of foam and plastic which forms the cam's shape to 
the design wavelength. 

Figure 8-2 shows a proposed solution for structural 
Bost reduction. 

The size of the cylindrical steel core is governed 
by volume requirements of conversion and generation apparatus. 
Let D = diameter of steel core 


core 


D = diameter of water surface cylinder 


skin 


The volume of the steel core per unit length = 


Z 


T e 
ie Doore 


Using the M.I.T. cam design, the volume of the 


shroud encompassing the steel core per unit length = 1.15% 


2 2 

Deskin ~ 74 Deore’ 

Let ey = conversion system volume density (Mm? /KW) 
Koen = generation system volume density (M>/KW) 

= 2 ; 
Kollass = S/M” of fiberglass coating 
_ 3 

Roly = S/M of polyurethane foam filler 


ey a = power to conversion system (KW/M) 





Polyurethane 
Foam Filler 


9 


skin /) 


Fiberglass 
surface 


Cylindrical Machinery 
steel Core 


Figure 8-2: Provosed Reduced Cost Structural 
Design. 
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The required core diameter is 


| > 


_ Nye 
Deore ~ 'F Peonv ‘conv * "conv Keen) ! 


Typical values of the above constants are, 


7 3 
SOnY = .23 M /KW 
f= .19 Mew 
gen : 
“conv = 92 


It was shown that the cost of a steel structure is 


199.71 D¢ 


Padel 


Therefore, 


Mt Pconv 
S mn lea ll 


core Pael 


mmne cost of the foam filler 


Z 


Cig S a 


| 


Kooly core 

















Lag) 


The cost of the fiberglass surface = 


KGlass Bode Tl sisi 
Pypically, 

K = 28 $/M> 

poly 

K = 4,3 $/M" 

glass . 


The total material cost of the shroud = 


2 - 22 p* «+ 15.5 D 


32.2 Dskin core skin 


Assuming that fabrication costs equal material costs (same 


as for all steel structure) 


_ — 2 
= 2(E25B ng. > 22 a. CCS Ie 


S ) 


shroud 


64.4 D + 31D 


I 


Zeer 


skin skin conv 


and the total cost of the structure is, 


2 
; _ 70 ae + 64.4 Desi + 31 vee sie 
structure Poel 

(9) 


- 
- : 
. i <a _ 
x 
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Revised versions of figures 6-5 and 6-15 may now be 
developed for optimizing structural costs using equation 9. 
Substituting parameters and non-dimensionalizing yields, 


70 (= \ + 64.4 (== y* 4 31(* 


) 
Say o. max max max 


3 Sr (10) 
max 65.4 














‘Equation 10 is plotted in figure 8-4 for the case 
with separation and in figure 8-3 for the case without 
separation. 

In figure 8-3 the minimum cost structure without 
separation occurs at yg eeeres = .66 with EEE ae = .79. 


‘In figure 8-4 the minimum cost structure with separa- 


flommoccurs at D/D,, -60 with P/P., =e ae 


Eo 7 ax 


It was shown in section 8.4 that the maximum output 
diameter as) occurs at .2034 Tr (without separation) and 
at .1383 T° (with separation). 


Total costs may now be expressed in terms of T, with 


the added 10% mooring cost. 


S struct (Without separation) = 
1.1[70 P +64.4[(.66) (.2034)]-“7* + 31(.66) (.2034) T°] 
i conv ; ° ° : : 
Pael 
4 2 
- eee tr. 28 Th + 4,58 7 
i <a aan 230 oa 


del 





BAe sane. ae 


6 





Zo 


02 4 06 AS 1.0 


D/D wax 


Figure 8-3: Minimum Cost Diameter Selection - No 


yeparation. 





ee And 2 ceo 


Laie 


$ /KW 


24 





D/D 


max 


Figure 8-4: Minimum Cost Diameter Selection - 


With Sevaration. 











ore 


=) 


(with separation) = 


Seat 
1.1{70 P + 64.4[(.60) ( 1383)]*r? + 31(.60) ( 1383) T°] 
. ey : : ; . 
ae 
77 Poo + 249 op? + 2.83 1° 
del 


8.9.1 Example fuel saver calculation with steel/foam 


construction. North Atlantic, with separation. 


8.5 sec 


+31 
I 


tol 
I 


37.1 KW/M 


O/T See cule 


eet 83 (figure 6-11) 


design = 38 


No = .32 
Pee tine = 11-9 KW /m 


Eee 2083 Po = eee 


Using equation 12, 


* Ste eonauenee ge $/KW 
_ 283 - 
MF ocr = 355 (PR) yields, PR = .52 


“band — -60 








oe 


@Ge— 377 + 369.(.52) = 579 $/KW 
AV-, = -85 NWUg = 51 
psy = 2224-2337) . 3.03 ¢/xWwh 


eT. os (ieee 


This compares with 4.49 ¢/KWh for the all-steel 


construction. 


for all six deep ocean locations around the U.S., 
steel and steel-foam construction, 


separation effect. 


Fuel saver analysis calculations were performed 
for all- 
with and without the 


The results are displayed’ in table 8-l. 


TABLE 8-1 


BREAKEVEN FUEL SAVING VALUES FOR VARIOUS WPG CONFIGURATIONS 


North 
Atlantic 
Mid 
Atlantic 


south 
Atlantic 


North 
Pacific 
Mid 
Pacific 
South 
Pacific 


Without Separation With Separation 
All-Steel jSteel-Foam | All-Steel | Steel-Foam 
T(sec) | P(KW/M)| nq |¢/KWh) ny ; ¢/KWhi] no | ¢/KWh| no | ¢/KWh 
J 
Sos Spell 20 lan. 5 25 | x25 | peels 5.0 
7.9 2526 ao 3. lela" | Ore lies | UBS iam Pec i Bel 
Gai 22:1 sls | oy ees | es | Oe 50). 31 | DEG 
Wi 8 81.0 22 a .29 , opamaieeeo 225”) 36 | sea 
| 
HG). 3 be. 0 21 a5. (Als ge eer Ono | BaD 
15) 42 2510 ah 7 532 | 27215. 32 | koe ea | LO 
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5.10 Analysis of Costs 

From table 8-1 the Middle Atlantic region and the North 
Pacific region yield comparable fuel saver values. The 
penalty for selecting high mean period (T) sites is signifi- 
cant: for a constant fuel saver value, period increases from 
7.9 to 11.0 seconds while power density increases from 22.1 
to 81 KW/M. 

The use of steel-foam instead of all-steel construction 
exhibits reduced costs, but a note of caution must be applied: 
Davis (ref. 7) showed that the cam's natural period of 
oscillation may be effectively raised by increasing the level, 
of inertia of the cam. ‘The steel-foam construction may be 
incapable of supporting the ballast at cam extremities for 
high levels of inertia. 

The separation effect provides the most dramatic 
reduction in costs. For most ocean sites this reduction in 
cost 1S approximately a factor of two. 

Using table 8-l and figure 8-1 for IHR = 11000 BTU/KWh, 
the lowest breakeven price of oil is 14 S$/BBL, corresponding 
to FSV = 2.5 ¢/KWh. It is unlikely that the WPG will ever 
achieve a breakeven cost lower than this. A median FSV = 
5.5 ¢/KWh corresponds to a breakeven oil price of about 


BO S$/BBL. 
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The WPG is a resource competing not only with conven- 
tional plants but also with other fuel-less generation 
systems like Wind Power and Ocean Thermal Energy Conversion. 
Costs of these competing resources have been estimated 
feets. 5, 6, 9, 13, 16, 21) to have FSV's ranging from 2.8 
to 3.5 ¢/KWh. For the WPG to compete in this range it must 
be configured as follows: 

1. Location in the open-ocean (approximately 100 
miles off-shore) with seas characterized by short 
mean wave periods and/or very high power densities; 

2. Separation between elements must be feasible; and 

3. A construction less costly than all-steel must be 
utilized, or a means of reducing the effective 


required diameter must be found. 
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CHAPTER 9 


CONCLUSIONS AND RECOMMENDATIONS 


The following list summarizes the most notable 
Semerusions reached in this revort: 
1. Ocean Waves 

a. Wave heights and periods are independent; 

b. Wave power distributions may be approximated 
by the exponential distribution; 

ce. East Coast waves are characterized by short 
wave periods and low power densities, while West Coast 
waves are characterized by long wave periods and high 
power densities; 

om The mean wave period remains essentially 
unchanged over the seasons; 

e. A typical ocean site disvlays avproximately 
mean annual power levels during the spring and fall, 60% 


of the mean in summer and 140% of the mean in winter. 


ee Modelling 

2. Simulation is a valuable tool for dealing 
with the complex WPG systen; 

b. For a given random ocean site there is one 
optimum cam diameter at a particular level of inertia; 


ce. Increasing the level of inertia reduces the 





a 


- 
- 


v 
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required cam diameter; 

d. the most cost-effective diameter is rated at 
68% of the mean power available without separation and at 
56% with sevaration for the all-steel structure. For the 
steel-foam structure the most cost-effective diameter is 
rated at 79% of the mean power available without separa 
tion and at 73% with separation; 

en “Lite parameter O,,/T may be used in determining 
the cam's hydrodynamic efficiency in a random sea; 

f. The separation effect results in higher effi- 
ciency, smaller maximum output diameter Dee) and 
); 


g. The most cost-effective level of energy storage 


smaller minimum cost/KW diameter ratio Oy 


for typical availabilities is 3 - 7 hours. 


3. Feasibility Studies 

a. wmeruenural “Costs are woronortional to tT, 

be. Direct mechanical conversion systems are 
infeasible; 

ce. Hydrostatic conversion is well-suited to the 
meG s 

d. DC electrical generation/transmission yields 
best results; 


e. The bandwidth of useable power may be ovutimized; 








iby, 


f. Flywheel energy storage is the most cost- 


effective storage mediun. 


4. Costs 

a. The annual fixed charge rate 1s approximately 
23% 3 

b. The concept of availability for the WPG is 
complicated by the variability of wave power; 

c. Steel-foam construction produces lower costs 
than all-steel construction; 

d. Breakeven costs with oil range from 14 to 


30 $/BBL. 


Recommendations 

Accurate open ocean data needs to be collected, 
Since power densities within 1 mile of the coast are too 
low to be economical. Data collection should concentrate 
on finding ocean areas characterized by nigh power densities 
and low mean annual wave periods. 

Model tests of free-floating multi-cam units should 
be conducted to determine (1) the feasibility of utilizing 
the separation effect, (2) mooring forces, (3) direction- 
keeping problems and (4) the limiting increase in natural 


period of oscillation with increases in levels of inertia. 
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Low cost structural concevts may be necessary to 
achieve competitive costs. Such designs should be con- 
ducted with prior knowledge of optimum levels of inertia. 

Inexpensive energy storage will increase the value 
of the WPG, but unit costs must be less than about 10 


$/KiWh to warrant storage use. 


The WPG will not likely prove to be the panacea 
of world energy problems. Under the most favorable 
assumptions it is barely competitive on today's energy 
market having a minimum breakeven oil price of 14 $/BBL;: 
under unfavorable assumptions it is infeasible with 
respect to today's energy costs. Before the WPG is elim- 
inated from feasible contention as an efficient energy 


producer, further investigation is required. 
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APPENDIX A 


MODEL DESCRIPTION 


FUNCTIONS 

1. NORML: Normal distribution for wave veriods 
(argument = RN). 

2. XPDIS: Exponential distribution for wave power 
(argument = RN). 

3. CAMEF: Cam efficiency (argument = 4D). 

4, GENEF: Generator efficiency (argument = actual power 
rating seen by generator). 

5. MULTI: Multiplier for increased power due to cam 


Separation (argument = cam separation). 


SAVEVALUES 

1. DEMND: Demand (KW). 

2. POWRT: Design generator power rating (KW). 
cee oie, Liz Not used. 

6. Pump efficiency (%). 

8. Mean power level (KW). 

9. Standard deviation of wave period (seconds). 
10. Mean wave period (seconds). 

12. Counter for assembling power from ten cams. 


13. Remaining storage capacity (KW). 


oo 
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14. Total storage cavacity (Ki). 
15. Accumulator of vower for delivery to generator. 


16. Counter for vower sources. 


VARIABLES 

1. Lamba: A/D (wavelength/diameter) 

2. POWER: Power available to each cam = power/length x 
length x MULTI. 

3. POWLN: Power/length. 

4, PERID: Wave period. 

5. PRATE: Power rating actually seen by generator. 

6. OUTGE: Power outage quantity 

7. PMPLS: Power available after conversion losses. 

8. CAMLS: Power available after cam losses. 

9. GENLS: Final power delivered by system. 

10. SUMPO: Power/length delivered by 10 cams. 

11. POWs Demand counter. 

Zee EMeC:s Counter. 

13. AVAIL: Power available in the storage. 


14. KAYD: Kd (wavenumber x separation distance). 


MATRIX SAVEVALUES 
1. Array (1x10) of cam lengths (m). 


Pine O) “of scam diameters (m). 
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Beeearray (1x10) of cam separations (m). 


TABLES 

fee eODs Period distribution. 

Ce POAVA: Power in wave distribution. 

3. LMDAD: “YD distribution. 

eee tos: Cam efficiency distribution. 

5. PULOS: Power leaving the conversion system. 
6. CALOS: Power leaving the cam. 

7. GELOS: Power out of the systen. 

8. OUTAG: Outage distribution. 


9. STOAR: Remaining capacity of storage. 
MODEL CONFIGURATION AND FLEXIBILITY 


The following model entities may be varied to test 
system response to different invut configurations: 
FUNCTIONS 
1. Any empirical cam efficiency distribution may be read 

into the model with CAMEF. 

2. Any empirical generator efficiency curve may be read 
into the model with GENE?F. 

3. Further model tests may yield different efficiency 
multiovliers than those currently in the model 


and these may be invout using the function MULTI. 
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MATRIX SAVEVALUES 

The model is designed to accent 10 different cam 
configurations for each trial run. The design length of 
the total system was arbitrarily assigned to be 1 Kil with 
each cam allotted 1/10 of the total system length. Cam 


diameters may be varied as desired. 


SAVEVALUES 

1. Demand: Power delivered and storage accumulation 
vary with changes in demand levels. 

2. Power Rating: The ability of the system to deliver 
power depends on the generator vower rating. This 
meper divided into the actual power available to 
the generator determines the input value for the 
generator efficiency curve. 

3. Conversion Efficiency: This is expressed as a constant 
and accounts for power losses after the cam but 
before the generator. 

4. Mean Power Level: Power levels vary from one location 
to another. This number along with the Function 
XPDIS comoletely describe the vower distribution. 

5. Period: The mean and standard deviation of wave veriod 
vary from one location to another. These variables 


affect optimum cam diameter, storage and nower 
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delivered. 
6. Storage Cavacity: The mean power delivered and the 
variability of this power are affected by changes 


in storage capacity. 


TABLES 
All significant model output is in the form of tables. 
Table means, standard deviations and cumulative vercen- 


tages are all imvortant identifiers of system output. 
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PENDIX B 
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